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A  calculated  procedure  will  be  developed  Sy  which  the  water  can  be  followed  as 
it  travels  downstream.  This  is  actual  water  routing  in  contrast  to  the 
classical  "water  routing  methods"  which  follow  the  hydrograph  wave  downstream. 
The  procedure  is  expanded  to  Include  the  tracing  of  conservative  water  chemis¬ 
try  as  the  water  moves  downstream.  Thus,  at  a  downstream  station,  it  is 
possible  to  predict  the  conservative  water  quality  parameters  and  compare  these 
parameters  with  measured  data  if  available.  Good  comparisons  between  calcu¬ 
lated  and  measured  conservative  substances  such  as  chlorides  and  conductivity 
indicate  the  calculat lonal  procedure  is  reasonably  accurate.  A  further  check 
on  accuracy  is  obtained  by  performing  a  mass  balance  on  the  total  water  volumes 
involved.  Predicted  concentrations  of  substances  such  as  total  phosphorus  and 
suspended  solids  do  not  compare  well  wlrh  the  measured  values,  indicating  that 
these  are  not  conservative  and  that  possible  deposition  and  resuspension  occur 
in  various  stream  reaches  during  the  storm.  Furthermore,  the  methodology  per¬ 
mits  the  estimation  of  the  fractions  of  various  parameters  at  a  given  down¬ 
stream  station  witich  have  been  derived  from  various  upstream  sources.  For 
example,  the  fraction  of  chloride  passing  the  Fremont  station  on  the  Sandusky 
River  which  was  derived  from  the  Honey  Creek  basin  can  be  calculated.  Thus,  if 
chloride  reduction  were  desired  for  Lake  Erie,  the  principal  upstream  sources 
can  be  identified  by  the  methodology  presented  herein.  - 


^TERIAL  TRANSPORT  IN  RIVER  SYSTEMS 
DURING  STORM  EVENTS  BY'VATER  ROUTING 


Water  QuaHty  Section 
Aray  Engineer  Dlatrlct,  Buffalo 


Frai^Verhoff 
Weet  VIrgTnIa  Unfverelty 
Morgantown,  West  Virginia 


fTis  eRJua 
DOC  TAB 
UoAnnounc  ed 
Justification 


Lake  Erie  Waatewater  Management  Study 
U.  S. I  Aray  Engineer  District,  Buffalo 
^  1776  Niagara  Street 

Buffalo.  NY  14207 


Avail  ai)d/or 
special 


TABLK  OF  CXMUTKMTS 


Siibj^pct  Page 

Lint  of  Fluure^  III 

Lint  of  Table*  vl 

Abstract  vll 

Introduction  I 

Cl  Iriilat  ional  Procedure  OevclopMent  1 

Water  Kihh  Ha  lance  1 

liydro)'rapb  and  Unneaaured  Inflow  Calculation  5 

Water  Routing  8 

Streaa  Chemlatry  Hass  Balance  9 

Flow  Fractions  ll 

Water  Chemistry  Concentration  11 

Water  Chemistry  Fractions  13 

Calculat tonal  Procedure  14 

Water  Transport  Calculations  17 

Concentration  Calculations  22 

Discussion  of  P.rrors  24 

Cl  Iciilat  tonal  Results  24 

Sandusky  River  Weir  Upper  Sandusky  33 

Sandusky  River  Near  Mexico  46 

Sandusky  River  Near  Fremont  46 

Llnearlxat Ion  of  the  Arca-Dlscharge 

Rt* lat  I onshl ps  46 

Time  Varying  Impacts  57 

Instream  Processing  of  Phosphorus  57 

Total  Water  Balance  Check  on  Accuracy  63 

Total  Phosphorus  Balance  63 

Management  Implications  67 

Conclusion  67 

Siimisary  67 

References  69 


ll 


FIGURES 


Number  Description  Fege 

1  Example  of  an  Area-Dlacharge  Curve  2 

2  A  Finite  Difference  Solution  Fora  for  Flow  4 

I  A  Llnearltatlon  of  Flow  or  Discharge  Area 

Between  Stations  6 

4  Flow  Fractions  10 

5  Chemistry  Fractions  12 

b  Example  of  Calculated  and  Approximated 

Unmeasured  Inflow  15 

7  Example  of  Hydrograph  Hatching  16 

8  Linear  Flow  Solution  Approximation  18 

9  Lower  Bound  on  Approximation  20 

10  Map  of  the  5>anduaky  River  Hatershed  23 

II  Example  of  Hydrograph  Mismatching  26 

12a  Flow  at  Upper  Sandusky  27 

I2h  Total  Phosphorus  st  Upper  Sandusky  28 

12c  Nitrite  and  Nitrate  at  Upper  Sandusky  29 

12d  Suspended  Solids  at  Upper  Sandusky  30 

12e  Chloride  at  Upper  Sandusky  31 

12f  Conductivity  at  Upper  Sandusky  32 

13a  Flow  at  Mexico  34 

13b  Total  Phosphorus  at  Mexico  35 

13c  Nitrite  and  Nitrate  at  Mexico  36 


111 


NiMib«r 


lid 

He 

nr 

14a 

I4h 

14c 

I4d 

14e 

I4f 

IS 

17a 
I7b 
I  Ha 
IRb 

18c 

19a 

19b 

1 9c 

20 

21 


View  of  One  fUiy  Flow  Using  Ltnesr  Ares-DlschArge 
Curve  Lj 

View  of  One  Osy  Flow  Using  Llnesr  Ares-DlschArge 
Curve  L2 

View  of  One  Dsy  Conductivity 


View  of  One  Day  Conductivity  Using  Linear-Area 
Discharge  Curve  L] 

View  of  One  Day  Conductivity  Using  Linear-Area 
Discharge  Curve  L2 

Tlee  Varying  Conductivity  Inflow  Concentration 

Conductivity  Results  Using  Variable  Inflow 
Concentration 


Page 

37 

38 

39 

40 

41 

42 

43 

44 

45 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

58 

59 


Iv 


D«scrtptlon 

Tlac  Varying  Total  Phoaphorua  Inflow 
Concentration 

Total  Phoaphorua  Raaulta  Ualng  Variable  Inflow 
Concentration 


Inatreaa  Proceaalng  of  Phoaphorua 


TABLES 


Nu»ber 

1 

2 

3 

4 


Docrlptlon  Page 


UngMgrd  Lateral  Inflov  Cheniatry  Concentrations 
for  tlie  Sandiiaky  River,  Storai  Beginning 
77/02/22  2S 

Conparlaon  of  the  Total  Water  Volunea  Proa 

Various  Sources,  Storm  Beginning  77/02/22  64 

Exaaple  of  a  Detailed  Total  Phosphorus  Budget, 

Storm  Beginning  77/02/22  65 

Example  of  a  General  Total  Phosphorus  Budget, 

Storm  Beginning  77/02/22  66 


vl 


ABSTRACT 


Flow  and  wnter  qiuillty  dat a  Ing  tn  197^hav»  bean  collerted,  analysed, 

and  plotted  at  nine  USCS  gaging  atationa  in  the  Sandusky  River  of  northern 
Ohio  for  aeveral  atora  eventa.  A  procedure  Itaa  been  developed  whereby  the 
water  as  it  travels  downstreaa  Is  followed.  The  calculations!  procedure  la 
further  expanded  to  Include  the  tracing  of  conservative  water  chealstry  as  It 
aoves  downatreaa.  Thus,  at  a  downstreaa  station  It  la  poaalblo  to  predict 
such  water  quality  paraaeters  and  coapare  thea  with  their  aeaaured  values. 

The  conservative  substances  such  as  chlorides  and  conductivity  coapared  well* 
Indicating  that  the  calculations  were  reasonably  accurate.  Such  substances 
as  total  phosphorus  and  suspended  solids  do  not  coapare  well  Indicating  these 
are  not  conservative  and  that  deposition  and  resuspension  occur  In  vsrlous 
reaches  of  the  atreaa.  Further,  the  aethodology  peralts  the  estlaatlon  of 
the  fractions  of  these  substances  at  a  given  downstreaa  station  which  are 
derived  froa  various  upstreaa  sources.  Various  factors  such  as  point 
sources,  river  flow  rate,  etc.,  are  discussed  as  to  how  they  affect  the 
transport  of  total  phosphorus,  suspended  solids,  chlorides,  conductivity,  the 
nitrogen  species,  and  other  water  quality  paraaeters. 


vll 


INTRODUCTION 


In  order  for  pl.«nners  .ind  i*nRtnecr«  to  he  .ihle  to  nn.tlyzo  .ind  evAliMte 
Instreo*  cnndltlonn  .«nd  recoiwend  water  quality  laanaxement  program*,  a 
qnoat  Inn  might  arlne  as  to  from  what  part  of  a  watershed  the  ntitjor  amount  of 
a  nutrient  or  flow  originates.  What  Is  needed  Is  a  simple  ralculatlonal  pro¬ 
cedure  requiring  a  minimal  amount  of  field  data.  The  object  of  this  report 
Is  to  present  such  a  procedure  developed  around  the  small  amount  of  data 
aval  I  able. 

\  c.i  Ico  lat  Iona  I  procedure  will  he  developed  by  wlilch  the  water  can  he 
followed  as  It  travels  downstream.  This  is  actual  water  routing  In  contrast 
to  the  classical  “water  routing  methods"  which  follow  the  hydrograph  wave 
downstream.  The  procedure  Is  expanded  to  Include  the  tracing  of  conservative 
water  chemistry  as  the  water  moves  downstream.  Thus,  at  a  downstream  sta¬ 
tion,  It  is  possible  to  predict  the  conservative  water  quality  parameters  and 
compare  these  parameters  with  measured  data  If  available.  Good  comparisons 
between  calculated  and  oM^asured  conservative  substances  such  as  chlorides  and 
conductivity  Indicate  the  calculations!  procedure  Is  reasonably  accurate.  A 
further  check  on  accuracy  Is  obtained  by  performing  a  mass  balance  on  the 
total  water  volumes  Involved.  Predicted  concentrations  of  substances  such  as 
total  phosphorus  and  suspended  solids  do  not  compare  well  with  the  measured 
values.  Indicating  that  these  are  not  conservative  and  that  possible  deposi¬ 
tion  and  resuspenslon  occur  In  various  stream  reaches  during  the  storm. 
Furthermore,  the  methodology  permits  the  estimation  of  the  fractions  of 
various  parameters  at  a  given  downstream  station  which  have  been  derived  from 
various  upstream  sources.  For  example,  the  fraction  of  chloride  passing  the 
Fremont  station  on  the  Sandusky  River  which  was  derived  from  the  Honey  Creek 
basin  can  be  calculated.  Thus  If  chloride  redjictlon  were  desired  for  Lake 
Frle  the  principle  upstream  sources  can  he  Identified  by  the  siethodology  pre¬ 
sented  herein. 

CALCULATIONAL  PROCEDURE  DEVELOPMENT 


Water  Ma^s  Balance 

Consider  the  mass  balance  for  stream  flow  Q  for  an  Infinitesimal  stream  ele¬ 
ment  X  between  large  point  source*  such  as  tributaries.  The  rate  of  change 
In  flow  for  the  element  I*  equal  to  the  inp<its  minus  outputs.  The  inputs 
consist  of  additions  from  small  point  and  diffuse  sources  along  x,  Che 
upstream  flow,  and  ground  water  Inflow.  Outputs  consist  of  stream  outflow 
from  the  element  plus  losses  to  the  stream  bank.  r.alns  or  losses  through 
ground  water  and  stream  banks  will  he  assumed  to  be  negligible.  In  reality, 
they  should  be  minor  during  storm  events.  Also,  they  are  extremely  hard  to 
approxlsMte.  The  differential  equation  takes  the  form 

-  q&a  ♦  Q  -  (Q  ♦  U  a,)  ...Eq(l) 


I 


where  A  •  strean  cross  secttonsl  flow  area 
Q  ■  discharge 

f  "  1 1  me 

X  ■  ilistanee  downs!  re. im 

rtx  •  ( nc remen t .« I  distance 

q  “  dlflnse  lateral  Inflows  jier  unit  length 

Point  sources  haslcally  represent  discontinuities  In  the  river  system  and  may 
be  Iwindled  In  two  different  ways  If  a  numerical  solution  technique  Is  used. 
Plrst,  point  sources  could  be  added  In  during  the  solution  procedure  by 
assuming  that  .tny  which  <vcur  over  an  Infinitesimal  reach  Ax  can  he  distri¬ 
buted  over  the  re.ich  or  secondly,  the  Ax's  could  he  chosen  such  that  the 
point  sources  .ire  added  to  the  stream  at  the  end  of  a  calciilat  lonal  step, 
meaning  the  flow  at  the  end  of  one  reach  Is  updated  for  the  beginning  of  the 
next  reach. 

Combining  terms  In  Kquatlon  1  and  taking  the  limit  as  Ax-»0  gives 


3A 


3t 


i!2  - 


3x 


,..F.q(2) 


This  Is  the  st.indard  differential  equation  used  .is  a  basis  for  all  water 
(hydrogr.iph)  routing.  Since  this  equation  has  two  unknown  parameters  Q  and 
A,  a  relationship  known  as  the  area-discharge  curve  Is  used  to  relate  the 
two.  In  other  Inst.inces  a  determlnlst  Ic  force  balance  Is  employed. 


Q  -  f(A) 


.F.q(3) 


Tills  relationship  Is  either  In  the  form  of  Rquatlon  3  as  determined  from 
actual  field  data  or  as  an  empirical  relationship  such  as 


oA 


.  Eq  (  4  ) 


where  and  o  must  he  approximated. 


Roth  Equation  3  and  4  have  been  used  with  Equation  2  for  water  routing. 
Verhoff  and  Melfl  (1978)  and  Verhoff,  Melfl,  and  Yakslch  (1979)  have  used 
field  data  while  Huang  (1978)  uses  an  emperical  relationship. 


Figure  I  Illustrates  a  typical  curve  for  discharge  versus  wetted  ares  for 
siihcrltlcal  flow.  This  curve  is  typical  for  most  streams.  The  average  water 
velocity  for  a  flow  Q  Is  Q/A  and  la  shown  on  the  diagram  as  the  slope  of  a 
chord  connecting  a  point  Q  on  the  curve  and  zero.  The  average  hydrograph 
velocity  la  the  slope  of  the  line  tangent  to  the  curve  at  point  Q.  As  can  be 
seen  from  the  curve,  the  hydrograph  velocity  Is  always  greater  than  the  water 
velocity.  As  the  water  moves  downstream,  the  hydrograph  will  catch  up  to  the 


c^rry  U  a  lilntJ^ncr  dovnst ream,  and  then  leave  It  behind.  If  a 
reference  point  la  attached  to  the  hydrograph,  the  water  would  be  observed  to 
laove  fr<>t«  left  to  right  through  the  hydrograph. 

Kquat  Iona  2  and  1  or  4  can  now  be  aolved  for  one  of  the  tw*>  paramatera 

lasinalng  the  dlffoae  lateral  lnfli>wa  jmt  unit  length,  q,  la  known  or  can  be 

reiKtinahly  approximated.  Depending  iip«>n  the  Information  wanted,  two  dif¬ 
ferent  approachea  mtty  be  taken.  If  q  la  either  known  or  can  he  reasonably 
approxI'Sitted,  Rquatlnna  2  and  1  can  he  solved  ualng  a  finite  difference 
scheme  giving  A  or  Q  as  a  function  of  x  and  t.  Such  a  solution  would  take 
the  form  In  Figure  2.  Figure  2  illustrates  a  solution  matrix  for  flow  for 
distance  downstreim  from  an  upstream  station  and  time  Into  a  storm  event. 

F.ach  nodal  point  which  la  Indicated  by  Q  represents  a  Solution  to  Equation  2 
for  equal  steps  In  tllstance  and  time.  Tributaries  that  occur  along  the  reach 
must  he  added  In  as  il  1 scont 1  nut t y  Jumps.  The  values  of  Q  at  time  xero  for 
all  iMstance  are  the  houniJary  conditions  while  Q'a  for  distance  xero  for  all 

t  lm«’  are  t*»e  Initial  conditions.  Values  other  than  those  calculated  for 

equil  distance  an.l  time  steps  must  h«*  Interpolated  from  these  solutions.  Tn 
ftils  case,  a  re lat  1  i>nsh I p  for  the  area-discharge  curve  would  he  needed  at 
each  X  distance  downstream.  For  «H  practical  purposes,  this  relationship 
could  h.*  approximated  from  the  relationship  at  end  points  of  the  stream  reach 
of  Interest. 


Hydri>graph  and  Unmeasuro«l _ ^nflow  Calculation 


If  lnform.it  Ion  Is  Just  neeiii'd  -it  the  .'nd  points,  an  alternative  approach 
could  he  taken  by  lntegr.it  Ing  F.quat  Ion  2  over  the  total  reach  distance 
g  I  V  I  nj; 


F.q(5) 


where  Lj  la  the  upstream  st.it  Ion  distance 

1.2  la  the  downstre.im  st.it  Ion  distance 


This  gives 


where  Q,  A, 


and  q  are  functions  of  time. 


...Eq(6) 


This  equation  Is  very  similar  to  the  "classic"  routing  techniques  which 
state:  the  change  In  storage  la  equal  to  the  Input  minus  output.  The 
Integr.il  of  are-i  Is  the  storage  between  stations  1  and  2  and  the  integral  on 
the  right-hand  aide  Is  the  nonpoint  Inflow  for  the  section.  The  stations  of 


S 


Intorrst  iir**  chosen  relatively  close  together  such  that  the  discharge  and/or 
flow  area  Is  approximately  linearly  dependent  upon  distance  between  the  sta¬ 
tions  during  storms  as  shown  In  Figure  3.  Here  the  amplitude  of  discharge 
(or  area)  and  variations  are  small  relative  to  the  wave  length  of  the 
tiydrograph  wave  form  passing  through  the  system.  As  a  guide  the  stations  are 
chosen  such  that  the  distance  between  them  Is  much  smaller  than  the 
hvdrograph  wave  length,  so  for  calculat lonal  purposes  the  water  surface  pro¬ 
file  (area)  Is  aasumtHl  to  be  relatively  linear  between  stations  as  the  wave 
progresses  downstream  as  shown  In  Figure  3  hy  the  heavy  straight  line.  This 
Is  onlv  tr«a>  If  th«Te  Is  no  abrupt  change  In  the  area-discharge  charac¬ 
teristics.  This  excludes  the  possible  Interruption  In  the  surface  profile 
due  to  large  point  source  Inflows  such  as  major  tributaries.  At  periods 
where  the  hydrograph  Is  starting  to  rise  from  some  base  flow  or  where  It 
peaks  relatively  sharply,  this  approximation  may  give  some  error. 

Assuming  a  linear  approximation  In  area  can  be  used,  Equation  6  gives 

Q2  -  Qj  —(A.  +  A,)  -  M  ...F,q(7) 

2  at  * 

where  M  replaces  the  nonpoint  Integral  and  L  Is  the  distance  between  stations 
(I.2~Lj).  At  this  point  the  term,  M,  could  also  Include  any  large  point 
sources  along  the  reach  such  as  tributaries.  Thus,  at  appropriate  times,  M 
replaces  all  the  point  and  diffuse  source  Inputs  between  the  two  stations 
with  a  reach  distance  L.  This  equation  may  be  used  to  calculated  H  between 
the  stations.  It  should  be  noted  that  F/^ustlon  7  only  has  to  he  applied  at 
the  two  end  points  of  the  reach.  This  equation  can  now  he  written  in  expli¬ 
cit  finite  difference  form  and  solved  for  the  local  Inflow  (Including  or 
excluding  point  sources)  between  two  stations  using  the  measured  hydrographs 
and  Kqiiat ton  3  or  A  at  the  two  end  stations.  Equation  3  may  be  easier  to  use 
than  Equation  A,  since  this  relationship  is  actually  obtained  from  field 
data.  This  Is  especially  true  since  hydrographs  are  needed  and  area- 
discharge  data  are  usually  obtained  In  the  process  of  defining  the 
hydrographs. 

It  must  b«'  remembered  that  M  Is  a  calculated  unmeasured  Inflow  dependent 
upon  Equations  7  and  3.  In  reality  M  might  not  exactly  Identify  with  the 
actual  unmeasured  flow  If  It  could  he  determined  Independent  of  any  routing. 

Equation  7  may  be  rearranged  to  define  a  dotmstream  flow  Q2  giving 


Q2  "  f (^2^  "  ^1  ”  ^  A2)  ♦  M  ...Eq(8) 

This  equation  may  be  used  to  "route"  a  hydrograph  from  an  upstream  station  to 
a  station  downstream  If  the  upstream  hydrograph  Is  known  along  with  the  area- 
discharge  curves  and  the  Inflow  function  M(t).  Thus  the  only  unkno«m  Is 
A2  which  can  he  found.  Q2  can  he  determined  from  A2  by  the  area-discharge 
relationship. 


A  nuMorlcal  technique  can  be  used  to  solve  Equations  7  and  8«  For  axaapla, 

In  explicit  finite  difference  froa  Equation  R  takes  the  fom 

"  ^.J  +  1  ^2,J 

■  '^2.j  ...Eq(9) 

where  1  denotes  the  upstream  station,  2  the  downstream  station,  and  j  time* 
Note  the  solution  Is  In  terms  of  area.  Areas  must  then  he  converted  to  flow 
throiif^h  Che  tise  of  the  area-discharge  curve. 

Q2  j  Is  the  discharge  at  the  >lownstream  station  at  time  J.  This  value  Is 
obtained  from  A2^j  using  the  area-discharge  relationships.  The  value  of 
A2  I  was  determined  by  the  aS'>ve  equation  for  the  previous  time  step. 

Inttlally  ^2 ,  J  ^2«J  *'*•*  ^  known  as  time  equals  zero. 

If  flow  and  consequently  the  flow  area  are  known  as  the  upstream  and 
downstream  station  and  the  ungaged  lateral  Inflow,  M,  Is  known.  Equation  9  can 
be  rearranged  to  solve  for  M  giving 


L 

2  it 


^2,J-H 


^2, j  ^  ^ 


...Eq(lO) 


A  quadratic  approximation  to  Che  area  as  a  function  of  distance  In  the 
Integral  which  was  linearized  above  has  been  attempted.  This  does  not  yield 
any  slm,>le  formula  or  calculat lonal  procedure  for  estimation  purposes. 


Water  kou t Ing 

The  tracing  of  the  water,  and  thus  Its  chemistry  as  It  moves  downstream,  can 
he  derived  from  Its  average  water  velocity.  If  the  velocity  Is  known  for  all 
times  diirlng  the  storm,  the  water's  position  la  known.  Thus,  for  a  given 
water  volume 


^  ■  Q/A  ■  V  ..,Eq(ll) 

dt  w 

where  Q  “  average  flow 

A  •  wetted  cross  sectional  flow  area 
a  •  axial  position  of  a  wmter  volume 

t  •  time 

v,,  •  average  water  velocity 


8 


If  In  aHBumed  constant  over  a  small  step  In  time  and  distance,  Equation  11 
In  Its  approximation  form  can  be  used  to  give  the  position  of  the  water  par¬ 
cel  as  It  travels  downstream.  Thus 


As  ■  v^At 


...Eq(l2) 


(>ncf  the  position  of  the  water  is  known  as  It  travels  downstream,  the  chem¬ 
istry  then  can  be  computed. 


Stream  Chemistry  Ma a s  Balance 


The  one  dimensional  longitudinal  mass  Iwilance  for  chemistry  for  a  differen¬ 
tial  element  gives 


at 


■f 


+  s  ♦ 

ource 


Ink 


...Eq(n) 


where  C  Indicates  chemistry  concentration 

D  Indicates  the  dispersion  coefficient 

Source  Indicates  chemistry  sources 
Indicates  chemistry  sinks 

The  above  equation  assumes  the  chemistry  Is  well  mixed  In  the  cross  section. 
Baker  (1979)  and  Baker  and  Kramer  (1979a,  1979b)  have  shown  good  mixing  to 
have  occurred  at  the  sampling  stations  used  on  the  Sandusky  River.  Generally 
the  smaller  the  stream  width,  the  better  the  mixing.  Also  there  is  better 
mixing  during  unsteady  storm  flo%n,  the  periods  for  %rhlch  the  analysis  was 
developed. 


By  traveling  with  the  water  as  It  moves  downstream  the  convective  term 

V  3c/3x  drop*  out. 

w 

The  first  term  on  the  right  hand  side  which  Is  the  longitudinal  dispersion 
term  is  assumed  negligible  for  the  following  reasons.  The  Influence  of 
dispersion  relative  to  convection  decreases  with  increasing  flow  rates. 
Although  dispersion  Is  of  much  less  Importance  durng  storm  events,  some 
researchers  even  neglect  It  during  steady  low  floM.  There  have  not  been  any 
dispersion  coefficient  studies  conducted  for  storm  flow  conditions.  This 
would  Indicate  It  may  not  be  possible  to  even  estimate  a  good  value  for  storm 
flow.  Also,  by  Including  dispersion,  the  calculations  are  greatly  com  pli¬ 
cated.  This  might  make  them  unwieldy  for  most  applications.  Thus  a  simple 
technique  which  does  not  Include  dispersion  could  have  wide  applicability  by 
a  variety  of  Government  agencies  and  other  study  organisations.  Finally  the 
consequences  of  excluding  the  dispersion  term  can  be  tested  by  comparing  the 
predicted  versus  measured  results  for  conservative  substances  such  as  chlo¬ 
rides,  and  conductivity.  Anticipating  the  results.  It  can  be  said  the  pre¬ 
dictions  without  dispersion  compare  well  with  the  measured  values.  Thus,  It 
appears  that  the  dispersion  term  can  be  neglected. 


FLOW  MASS  BALANCE 


BrcrtuAt*  of  the  annuapt  Ions  made,  this  procedure  might  not  yield  "accurate*' 
results  under  conditions  where  dispersion  Is  Important*  The  "accuracy**  la  to 
be  defined  hy  the  user  from  the  results  for  conservative  substances. 

Thus  Filiation  11  becomes 


^  -  S  *  S 

at  °‘**‘C*  ...F,q(14) 

This  is  essentially  the  I.agranglan  form  of  the  mass  balance  on  the  substance 
i>f  Interest. 

This  equation  basically  says  that  the  substance  concentration  In  a  fixed 
volume  of  water  will  change  because  of  Internal  generation,  ^ource 
nitrite  generation  by  ammonia  oxidation),  or  because  of  internal  disap¬ 
pearance,  (e.g.  sediment  settling  to  the  bottom).  Since  this  analysis 

wan  based  on  conservative  substances  both  and  are  r.ero.  Thus  so 

long  as  there  are  no  Internal  lnp<its  for  mixing  the  concentration  remains 
constant  as  the  water  moves  <liiwrist ream.  Since  all  external  sources  are 
assumed  to  occur  at  a  point,  no  accumulation  is  possible.  Only  mixing  need 
b«*  considered  at  the  point. 

4s  in  example,  consider  a  steady  state  flow  of  K  with  concentration  B.  If 
another  sfiurce  Is  added  with  a  flow  of  C  and  concentration  D,  the  total  flow 
Is  A+C  with  a  resultant  concentration  of  ( AB+CD)/(A+C).  The  fraction  of  the 
mainstream  flow  Is  A/(A+C)  and  that  of  the  added  source  is  C/(A+C).  The 
fraction  of  the  chemistry  ( AB+CD)/( A+C)  that  originates  from  the  mainstream 
flow  Is  AB/(ABK:n),  while  that  from  the  added  source  Is  a)/(AB4CD). 

Flow  Fr ac t Ions 

Following  this  same  prnce<lure,  the  flow  fractions  at  time  1+1  for  the  gaged 
and  ungaged  flow  can  he  calculated  as  shown  In  Figure  4. 


T 


OC 


1+1 


Qj  qAs  +  Qj. 


+ 


Q 


c 


Qj  +  qAs  +  Qp 


...Faj(15) 


1  1 


Qj  qAs  4  Qj. 


qAs 


4  qAs  +  Qg 


...Eq(16) 
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CHEMISTRY  MASS  BALANCE 
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Fluure 


where  F  Indicates  fraction 
C  Indicates  gaged 
U  Indicates  ungaged  (flow  is 

distributed  as  lateral  Inflow) 
QO  Indicates  gaged  flow 
Qll  indicates  ungaged  flow 


Also  the  fraction  of  the  flow  from  the  upstream  Ith  element  at  the  l+l£h 
Pitlnt  is 

QC  Ql' 

1  -  F  -  ,  • 

i-*-l  i  +  1 


Water  Chemistry  Concentration 


By  doing  a  steady  state  mass  balance  on  the  chemistry,  the  chemistry  con 
cent  ration  for  the  water  at  time  l+l  can  be  shown  to  be 


C 


!♦! 


*  Qj. 


Ql  ♦  qAa  4 


C^.  qAs 

Ql  ♦  qAa  ♦  Op 


...Kq(l7) 


where  C  indicates  chemistry  concentrations 

On  the  right  hand  side,  the  first  term  represents  the  dilution  of  the 
upstream  chemistry  by  the  flow  added  from  point  (tributary)  and  nonpoint 
(ungaged  lateral  Inflow)  sources;  the  other  two  terms  represent  the  addition 
by  a  tributary  and  unmeasured  Inflow  for  the  Incremental  reach  %,  respec¬ 
tively.  Knowing  sufficient  Information  on  all  the  tributaries  and  each 
ungaged  inflow  area,  the  tributary  and  ungaged  concentrations  can  be  traced 
through  the  system. 

Water  Cheml stry  Fractions 

The  chemistry  due  to  various  inputs  Is  not  calculated  directly,  but  Is  traced 
as  a  fraction  of  the  concentration  given  In  Equation  17.  As  shown  In  Figure 
6  the  gaged  and  ungaged  chemistry  fractions  are  given  by 


rCC  r 


T.CC  ^ 


Q,c^  *  C^qA.  4 


...Eq(18) 
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Q' 

Q  C 
i  i  1 


C^.qA^ 

QjC,  .  C|j,A.  ♦ 


...Eq(l9) 


where  CC  denotes  Rnxed  chemistry 

CU  denotes  unRSRed  chemistry 


Again  the  fraction  the  upstream  station  Is  1  -  F 


CC 

1+] 


-  F 


a‘ 

i-fi 


The  first  term  on  the  right-hand  side  represents  the  dilution  of  the  chem¬ 
istry  fraction  as  the  water  travels  downstream,  while  the  last  term  Is  the 
fractional  addition  due  to  an  ungaged  source.  Equations  15  and  18  may  he 
repeated  for  several  tributaries  and  other  gaged  point  sources  to  be  moni¬ 
tored.  Tributaries  which  are  not  monitored  are  grouped  with  the  nonpoint 
ungaged  sources.  These  ungaged  sources  may  be  broken  Into  as  many  aubsources 
as  desired  by  using  Equations  lf>  and  19  as  long  as  appropriate  values  for  C 
and  q  for  each  subset  exist  or  can  be  reasonably  approximated. 

Water  to  be  traced  starts  out  at  the  upstream  station.  Here  It  has  the 
measured  chemistry  concentrations  and  flow  at  that  station;  and  the  fraction 
from  the  upstream  station  Is  one  while  all  other  fractions  are  zero.  As  this 
water  travels  downstream,  the  Influence  of  the  upstream  station  becomes  lens 
and  less  while  the  Influence  of  the  unmeasured  contribution  Increases.  When 
this  water  hits  a  tributary  during  a  routing  Increment  x,  the  tributary 
Input  affects  the  water  at  this  point.  Thereafter,  the  tributary  fraction  Is 
dl  luted. 


The  calculations  translate  the  chemistry  and  water  to  the  downstream  station 
where  comparisons  are  made  with  measured  parameters.  The  preceding  deriva¬ 
tions  were  made  only  for  this  comparison  and  they  will  not  necessarily  give 
accurate  Information  bet%#een  the  stations  of  Interest. 


CALCULATIONAL  PROCEDURE 


The  fundamental  principles  useful  for  understanding  the  transport  of 
fliaterlals  In  streams  have  been  discussed  In  the  previous  sections.  The  goal 
of  this  section  Is  to  present  the  sequence  of  calculations  required  to  follow 
the  nutrients  and  other  materials  as  they  move  through  the  main  stem  of  a 
river  system.  Raslcally  these  calculations  are  accomplished  In  two  steps. 


lA 


SANDUSKY  RIVER  FROM 
BUCYRU5  TO  UPPER  SANDUSKY 
STORM  beginning  77/02/2Z 

—  calculated  local  inflow 

—  approximated  local  inflow 

INCLUDING  tributaries 


Figure  6 


SfiNDUSKY  RIVER  NEAR  UPPER  SANDUSKY 
STORM  BEGINNING  77/02/22 
—  OBSERVED  HYOROGRAPH 
—  calculated  HYOROGRAPH 
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First  the  n»<)s^ein<‘nt  nf  the  <#.iter  fron  the  upstreAm  AtAt  ton  And  the  tributaries 
to  the  downstreAn  stAtlon  must  be  traced.  Once  this  nK)vement  of  the  water  la 
quantified  the  transport  of  nutrients  can  he  calculated. 

Ttu*  calculations  described  Iwreln  trace  the  water  and  substances  as  they  are 
t  ransp»>rte<l  f  ri>m  an  upstream  station  of  a  river  to  a  downstream  station.  In 
order  to  trace  the  water,  some  relationship  between  the  discharge  and  the 
river  cross-sectional  area  (area-discharge  relationship)  must  he  used.  For 
exact  calculations  this  discharge  relationship  should  be  measured  as  a  fiinc- 
tli>n  of  distance  di>wnst  ream.  However,  for  calculat  tonal  purposes  l»eretr.  It 
will  h«“  assumed  that  the  area-illscharge  relationship  at  each  potul  is  an 
aver.ige  of  that  .it  the  upstre.im  .ind  downstream  stations.  Also,  for  calcula¬ 
tions  In  the  main  stem  of  a  river  It  Is  necessary  to  know  for  any  tributary 
of  Interest  Its  flow  and  material  Inputs  at  Its  mouth.  However,  the 
measiirliig  stations  on  the  tributaries  are  usually  located  some  distance 
upstream.  Tlierefore  the  tributary  Inputs  must  he  translated  downstream  to 
the  mouth  before  the  nwiln  stem  calculations  can  he  completed. 


l^ater_  Transport  Calculations 

Many  procedures  have  been  developed  for  the  solution  of  Equation  2  including 
finite  difference  schemes  and  the  method  of  characteristics.  The  suijor  goal 
of  these  solution  techniques  has  been  to  model  and  predict  the  height  of  the 
water  (hydrograph  wave)  during  a  storm.  Tlie  goal  of  the  paper  presented  here 
Is  quite  illfferent  and  thus  the  manner  In  which  this  equation  Is  used  will  he 
different.  ''Ince  the  water  flow  Is  known  both  at  the  upstream  and  the 
ilownstre.im  stations,  the  main  use  of  Fx^uatlon  2  Is  for  the  calculation  of  the 
unme.isure.l  l.iteral  Inflow,  q.  The  equation  Is  then  used  to  predict  the 
downst  rivi'S  hydrugr  iph  with  smixithed  values  of  the  calcul.ited  Inflow  q. 

The  equation  us«-d  to  cilculate  M,  the  local  Inflow  betw«>en  the  two  main  stem 
stations  Including  tributaries.  Is  given  In  F.quatlon  10,  where  the  length 
between  st.it  Ions  Is  T.  and  the  time  Interval  Is  ^jt.  Since  the  discharge  Q 
and  the  cross-sectional  area  A  are  known  as  a  function  of  time  at  both 
upstre.im  and  downstre.sm  stations  (the  area-discharge  curve  at  each  station  Is 
employed),  the  only  unknown  In  this  equation  for  any  time  Interval  Is  Mj+|. 
Thus  the  lateral  Inflow  for  that  section  of  the  river  can  he  calculated  as  a 
function  of  time  for  the  Increments  of  time.  Such  a  calculation  between 
Hucyrus  and  Upper  Sandusky,  Ohio,  on  the  Sandusky  River  Is  shown  In  Figure  b. 
The  Inflow  fiincton  Is  somewhat  Irregular  and  a  smoothed  version  of  this  curve 
Is  ilso  shown  In  this  figure. 

The  smiMitlieil  Inflow  function  M  obtained  from  Figure  6  la  used  along  with  the 
upstri’im  discharge  and  discharge  area  for  a  storm  along  with  the  past 
downstre.im  dlnclvirge  to  predict  the  future  downstream  discharge  area  and 
dlsclvirge  using  F.quatlon  9.  The  resulting  predictions  are  shown  in  Figure  7. 
It  can  be  si>en  that  using  the  smoothed  input  functions  causes  little 
deviations  In  the  downstream  hydrograph. 
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The  CAlciilHteil  Inflow  function  M  contalnn  both  the  meAHured  tributary  Inflow 
and  that  rexiiltlnK  from  unmeaniired  areax  draining;  directly  Into  the  main  atem 
of  the  river.  T'je  unmeaaureil  Inflow  la  calculated  hy  aubtractlng  the  tribu¬ 
tary  Inflow  from  the  calculated  total  Inflow.  Occaalonally ,  aome  calculated 
valiiea  of  the  net  <in«eaaurtfd  Inflow  are  negative.  Thla  could  he  reallatlc 
becaiiae  there  can  he  a  algnlf leant  recharge  of  the  groundwater  eapeclally 
durlitg  the  Hummer  i%)ittha  In  thla  area  (western  Ohio).  Howver,  since  this 
occurs  rarely  and  randomly  during  the  calculat  luna I  procedure,  the  unmeasured 
Inflow  was  constrained  to  he  equal  to  i>r  greater  than  lero. 

The  tributary  Inflow  which  Is  subtracted  from  the  total  Inflow  Is  that  eval¬ 
uated  at  the  tributary  mouth.  Since  the  mt'asurement  jvjlnt  on  the  tributary 
Is  some  distance  upstream  from  the  mouth,  It  Is  necessary  to  calculate  the 
inflow  at  the  mouth.  This  Is  accomplished  hy  using  Equation  10  to  calculate 
an  M  for  the  area  upstream  of  the  measurement  point.  The  upstream  station  Is 
taVen  as  the  upp«>rmoHt  headwaters  of  the  stream  with  a  rero  area-dlscharg.. 
Curve  and  hydrograph.  Tltls  local  Inflow  Is  then  prorated  hy  the  area  drained 
In  the  short  distance  from  tributary  measurement  station  to  Its  mouth  and 
tised  In  F.quatlon  *1  to  calculate  the  hydrograph  at  the  tributary  mouth. 

In  all  the  cases  studied  In  the  Sandusky  River,  the  Inflow  from  the  measure¬ 
ment  point  to  the  tributary  mouth  was  minor.  The  major  change  that  occurred 
In  translating  the  tributary  flow  from  the  measurement  point  to  the  mouth  was 
the  delay  In  time  at  which  the  flow  entered  the  main  stem  of  the  river. 

•\t  this  p«)lnt  the  Inflows  to  the  main  stem  of  the  river  are  defined  for  the 
c'otlre  storm  event  and  the  water  can  he  traced  as  It  progresses  downstream. 
The  basic  equation  used  to  specify  the  position  of  the  water  Is  Equation  11 
or  Its  approximation  F.quatlon  12.  Using  a  linear  approximation  for  the  flow 
between  stat  Ions  and  using  the  known  fli>ws  at  the  upstream  and  domstrearn 
Htal Ions  as  well  as  the  known  tributary  Inflows  it  Is  possible  to  predict  the 
dtsclurge  In  the  main  stem  of  the  river  as  a  function  of  distance  for  any 
Inst.inf  In  time.  This  Is  shown  In  Figure  R  for  a  time  tj  and  the  next  time 
Increment  t  -f^t  which  equals  tj+j.  Given  that  the  parcel  Is  at  a  point 
sj  at  time  tj  with  flow  rate  Qj,  w*’  wish  to  calculate  the  position  of  the 
parcel  at  time  tj+|.  From  the  area-discharge  relationships  at  the  end  sta¬ 
tions  the  cross-sectional  area  of  the  river  at  the  point  of  Interest  Is  esti¬ 
mated  and  the  water  velocity  Is  calculated  as  Q^/Ap  During  the  time 
Interval  t  •  ^l+l  “  ^i  Ibe  water  will  move  to  a  new  position  Sj^j.  This 
position  Is  calculated  by  the  following  formula. 


B 


l+l 


*  (Qj/A^)At 


.  .  .  Eq  (20) 


This  calculation  then  positions  the  water  parcel  at  distance  sj4.(  and  tj4.|« 
This  point  Is  Indicated  by  Qj+l  on  Figure  8.  This  process  Is  then  repeated 
with  the  flow  as  a  function  of  distance  drawn  for  time  t|4.2. 


19 


UISTANCJ  DOWNSTREAM 


For  Kreater  accuracy,  the  following  fornula  was  iiaed: 


*141  -  *1  ♦  4  Qj^,/Aj^j)/2}At  ...Eq.(21) 

Thla  equation  la  uaed  Iteratively  after  a  first  estimate  of  s^^.]  and  Q^-fi  la 
obtained  from  Figure  8  until  the  approximation  for  a^i  and  match  with 

the  flow  curve  at  t(4.i>  This  Is  shotm  In  Figure  8. 

This  calculation  process  Is  repeated  In  the  same  manner  except  when  a  tribu¬ 
tary  Is  crossed.  If  the  calculation  of  the  new  position  of  the  water,  *14.], 
Indicates  that  a  tributary  has  been  passed  the  calculation  Is  modified  such 
that  the  Increment  of  time  At  Is  divided  Into  two  parts.  The  flrt  part 
carries  the  water  to  the  tributary  using  the  velocity  calculated  at  the 
starting  point.  For  the  remainder  of  /It,  a  distance  of  travel  la  calculated 
using  the  main  stem  stream  flow  Just  downstream  of  the  tributary  mouth  as  the 
starting  point. 

This  calculation  then  permits  the  position  of  the  water  to  be  known  as  a 
function  of  time.  Usually  the  water  Is  started  at  the  upstream  station,  s  ■■ 
0,  at  different  times.  Some  water  Is  "tagged"  as  It  starts  downstream  before 
the  storm  begins  at  the  upstream  station.  Thereafter,  water  from  different 
parts  of  the  upstream  hydrograph  are  chosen  at  time  intervals  such  that  Che 
fate  of  the  upstream  water  can  be  followed  through  the  main  stem  of  the 
stream.  Thus  It  Is  now  possible  to  predict  when  water  from  a  poaltlon  In  the 
upstream  hydrograph  reaches  the  downstream  station.  Usually,  the  water 
travels  slower  than  the  storm  wave  (hydrograph)  and  hence,  the  water  at  the 
peak  of  the  hydrograph  at  the  upstream  station  reaches  the  downstream  station 
after  the  hydrograph  has  peaked  at  this  station. 

In  addition  to  computing  Che  poaltlon  of  the  water  It  Is  possible  to  deter¬ 
mine  the  composition  of  Che  water  which  arrives  at  the  downstream  station. 
That  Is,  the  fraction  of  the  water  which  resulted  from  the  upstream  station, 
the  tributaries,  and  the  unmeasured  Inputs  can  be  found.  The  equations  for 
these  calculations  were  given  as  Equations  15  and  16.  The  fraction  for  each 
of  the  Inputs  Is  updated  during  each  time  Increment,  l.e.  as  tlsM  goes  from 
t^  to  t(>|.  Thus  It  would  be  possible  to  calculate  the  fractlona  of  water 
from  each  source  as  a  function  of  distance  for  the  water. 

Occasionally,  during  the  calculations  the  discharge  predicted  from  the  linear 
curve  between  the  t««o  main  stem  stations  Including  discontinuities  caused  by 
the  tributary  Input  will  drop  below  that  of  the  low  flow  discharge  in  the 
river.  This  Is  shown  in  Figure  9.  In  this  Instance  the  flow  Is  smoothed 
along  the  low  flow  line  as  shown  In  Figure  9.  This  occurs  rarely,  but  when 
It  docs  occur  It  Is  during  a  rapidly  rising  or  declining  flow  of  a  storm. 

The  rapidly  declining  stage  of  the  storm  Is  shown  In  Figure  9.  It  should  be 
noted  that  this  sproxlmatlon  does  not  Introduce  much  more  error  since  the 
linear  approximation  for  the  discharge  as  a  function  of  distance  la  good 
during  these  rapid  changes  In  flow. 


21 


Coil ce ntratlon  Calculatlon« 

Thr  concentration  of  the  water  at  the  tipatream  station  Is  known  from  the  data 
measured  there.  As  the  water  leaves  the  upstream  station  It  mixes  with  water 
from  the  unmeasured  Inflow  areas  of  the  hasln  and  with  the  water  from  the 
tributaries.  Since  the  concentration  of  materials  In  the  unmeasured  wmter  or 
the  tributaries  Is  not  necessarily  the  same  as  that  of  the  water  In  the  main 
stem  of  the  stream,  the  concentration  of  material  In  the  water  changes  as  It 
proceeds  downstream.  Stated  simply,  the  water  In  the  main  stem  Is  diluted  by 
water  from  tributaries  and  unmeasured  areas  of  the  basin. 

If  the  storage  capacity  of  the  river  at  the  point  of  external  water  entry  Is 
small  compared  to  the  flow  rates  of  the  water,  the  steady  state  mixing  con¬ 
cept  can  be  used  to  calculate  the  resulting  concentration  In  the  water.  This 
assumption  Is  very  good  for  all  streams  which  do  not  have  large  natural  or 
man-made  Impoundments;  this  la  the  case  for  the  Sandusky  River  which  Is 
studied  herein,  "the  steady  state  mixing  equation  was  given  as  Equation  17. 

This  equation  Is  calculated  for  each  Increment  in  time  and  thus  the  con¬ 
centration  of  the  material  of  Interest  Is  determined  as  a  function  of 
distance  since  distance  Is  also  a  function  of  time.  It  must  be  remembered 
that  this  calculation  presumes  the  substance  Is  conserved  in  the  water. 

This  calculation.  In  particular,  yields  the  predicted  concentration  at  the 
downstream  measurement  station  on  the  main  stem.  At  this  point  the  calcu¬ 
lated  concent  rat  Ion  can  he  compared  with  the  measured  value.  It  la  expected 
that  conserved  substances  will  compare  favorably  and  that  non-conserved 
substances  will  exhibit  deviations  from  the  predicted  values.  These 
deviations  will  give  Indications  as  to  transport  processes  e.g.  If  the 
measured  concentration  Is  greater  than  that  computed,  there  might  have  been 
some  net  resuapenslon  of  that  substance  for  the  stream  bottom  and  banks  into 
the  flowing  water.  Thus,  the  comparison  between  calculated  and  measured  con¬ 
centration  at  the  downstream  station  is  valuable  In  that  it  permits  a  test  of 
the  accuracy  of  the  computations  for  conservative  substances  such  as  chloride 
or  conductivity  and  It  yields  Information  as  to  the  transport  processes  of 
nonconservative  substances  such  as  suspended  sediment  and  total  phosphorus. 

In  addition  to  the  cmputatlons  on  the  concentration  of  the  water  as  a  func¬ 
tion  of  distance,  It  Is  possible  to  calculate  the  fraction  of  the  substance 
In  the  water  which  comes  from  different  sources.  These  fractions,  from  the 
upstream  station,  from  tr ihvitar les ,  and  from  unmeasured  areas  are  calculated 
as  the  water  proceeds  downstream.  The  equations  for  these  calculations  were 
given  as  Equations  18  and  I9. 

It  must  be  remembered  that  these  calculations  are  for  conservative  subatancea 
such  as  chlorides  or  conductivity.  However,  the  fractions  of  other  sub¬ 
stances  can  contain  significant  Inforsuitlon  as  to  the  processing  of  these 
substances  In  the  stream.  For  example,  one  might  find  a  correlation  between 
decreased  quantities  of  a  given  substance  e.g.  total  phosphorus  and  the  time 
when  a  given  fraction  of  total  phosphorus  was  predicted  to  be  from  a  given 
source,  e.g.  a  particular  tributary.  Thun  one  might  conclude  that  the  total 
phosphorus  from  that  source  might  he  settling  out  In  the  main  stem  of  the 
river. 
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The  computational  technique  presented  herein  for  the  tracing  of  substances  In 
a  river  contains  a  number  of  expedient  assumptions.  These  assumptions  are 
made  with  the  recognition  that  the  quantlatlve  conclusions  derived  may  have 
some  error  associated  with  them.  However,  these  assumptions  are  expedient 
l.e.  the  computations  could  not  be  accomplished  without  them,  and  the 
understanding  of  river  transport  to  be  gained  from  calculations  outweighed 
the  errors  Incurred. 

The  N*«t  teat  of  the  compiitat  Ion  technique  Involves  the  comparison  of  the 
predictions  and  measurements  of  the  water  flow  and  the  conservative  substance 
concentrations  at  the  downstream  stations.  If  the  predictions  and  the 
measured  quantities  at  the  downstream  station  match  reasonably  well  It  will 
be  concluded  that  the  assumptions  made  do  not  adversely  effect  the  com¬ 
putations.  As  will  he  shown  later.  It  appears  that  the  predictions  and 
measurements  compare  well. 

It  Is  worthwhile  to  summarize  these  expedient  assumptions.  The  water 
discharge  Is  presumed  to  be  linear  as  a  function  of  distance  between  the  two 
main  stem  river  stations.  To  presume  any  higher  order  function  would  be  too 
complex  mathematically.  The  tributary  discharge  and  concentrations  were 
approximately  translated  downstream  from  the  measuring  station  to  the  mouth 
where  the  water  entered  the  main  stem  of  the  river.  It  Is  generally 
Impossible  to  measure  the  water  flow  at  the  mouth.  The  method  used  to 
translate  the  water  downstream  (basically  a  time  delay  and  area  flow 
augmentation)  were  about  as  simple  as  could  be  conceived.  Other  more 
sophisticated  calculations  could  be  used  for  this  purpose  depending  upon  the 
application  l.e.  the  distance  bet%tecn  the  measuring  station  and  the  mouth  of 
the  tributary. 

It  Is  also  assumed  that  the  unmeasured  Inflow  cannot  become  negative.  This 
assumption  Is  Invoked  Infrequently.  However,  It  may  not  necessarily  be  true 
since  a  negative  Inflow  from  unateasiired  areas  Is  possible  when  water  from  the 
river  enters  the  soil  for  ground  water  recharge.  Also,  the  flow  In  the 
stream  is  never  permitted  to  drop  below  that  of  low  flow  conditions.  Again 
this  condition  Is  Involved  Infrequently.  A  higher  degree  approximation  of 
the  discharge  as  a  function  of  distance  might  alleviate  this  problem. 

CALCULATIONAL  RESULTS 


The  calculations!  procedure  developed  above  was  used  to  study  storms  on  the 
Sandusky  River  basin.  Water  routing  was  started  at  the  USCS  gaging  station 
near  Bucyrus  with  the  calculated  results  being  compared  to  measured  data  at 
downstream  stations  near  Upper  Sandusky,  Nexlco,  and  Fremont.  A  map  of  the 
watershed  Is  shown  In  Figure  10. 

These  data  «iere  collected  and  analyzed  by  Heidelberg  College,  River  Studies 
Laboratory,  for  the  U.  S.  Army  Corps  of  Engineers.  See  U.S.  Army  Corps  of 
Engineers,  Buffalo  District  (1978a,  1978b,  1978c).  The  arsa-dlschargs  rela¬ 
tionships  were  developed  from  actual  field  data  collected  by  the  USCS  at  the 


Table  1  -  Ungaged  Lateral  Inflow  Concentrations 
Stora  of  77/02/22 


Ungaged  Lateral  Inflow 
Concentrations  Excluding  Trlbntarlea 


Station 

TP 

N02+N03 

«g/l 

SS 

m/i 

Cl 

■«/l 

Cond. 

uahos 

Bucyrua 

Upper  Sandusky 

0.2 

1.0 

40 

25 

650 

Hex  Ico 

0.2 

3.0 

75 

50 

500 

Freaont 

0.2 

2.0 

40 

50 

500 
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tlae  they  defined  or  redefined  the  discharge  relationships  at  the  stations 
used.  The  following  results  and  plots  are  just  to  give  an  Indication  of  what 
type  of  Inforaatlon  can  be  derived  frosi  this  procedure.  Figure  6  Illustrates 
a  calculated  and  approx laiated  ungaged  Inflow  Including  the  tributary  of 
Interest,  Broken  Sword  Creek.  This  H  Is  for  the  reach  of  Bucyrua  to  Upper 
Sandusky  for  the  store  of  77/02/22.  Figure  7  shows  an  excellent  aatch  on  Che 
downstreae  hydrograph  for  this  store. 

For  this  reach,  It  was  found  that  the  ungaged  Inflow  went  negative  at  soee 
periods  of  tlee  when  the  tributary  hydrograph  for  Broken  Strard  Creek  was 
subtracted.  Since  this  adjusted  calculated  M  and  the  hydrograph  for  Broken 
Sword  Creek  are  used  in  the  water  routing  to  define  the  cheelstry,  they  are 
added  together  to  define  a  new  N.  Figure  II  shows  the  resulting  alsMtch  In 
the  hydrograph  as  a  result  of  this  adjustaent. 

Figures  12a-t2f,  Figures  13a-13f,  and  Figures  14a-14f  show  the  results  for 
the  Upper  Sandusky,  Mexico,  and  Freaont  stations  for  the  stom  of  77/02/22. 
The  aatertal  concentrations  froa  unaeasured  areas  are  asauaed  to  be  constant 
over  the  tlae  period  and  are  given  In  Table  1. 

Sandusky  River  Wear  Upper  Sandusky 

Figure  12a  shows  the  unaeasured  Inflow  to  be  predoalnant  only  during  the  22nd 
and  23rd.  It  can  be  expected  to  have  relative  Influence  on  the  chealatry 
only  during  this  period. 

Figure  12b  shows  reauspenslon  of  total  phosphorus  throughout  the  duration  of 
the  atora  for  an  unaeasured  lateral  inflow  concentration  of  0.2  ag/1. 
Increasing  the  locsl  Inflow  concentrst Ion  to  a  high  value  would  only  have  a 
aa Jor  effect  on  the  first  two  days.  Tills  aay  reduce  the  resuspenslon  of 
total  phosphorus  during  these  days.  The  peak  that  shows  up  over  the  lot  and 
2nd  results  froa  a  spike  at  the  Bucyruo  station.  This  spike  Is  represented 
by  one  high  nisiber  In  the  data  for  this  station.  Since  one  aeasureaenC  Is 
Involved,  this  aay  Indicate  an  erroneoua  data  point.  However,  since  this 
saae  trend  shows  up  for  suspended  solids  and  the  station  Is  just  downstreae 
of  a  sewage  treataent  plant,  this  point  Is  probably  real. 

Figure  12c  generally  shows  a  loss  of  nltrate/nltrlte  nitrogen  froa  the  systea 
for  the  unaeasured  Inflow  concentration  which  was  aosuaed  In  the  coaputatlon. 

Figure  12d  shows  a  general  resuspenslon  of  sedlacnt  froa  the  bottoa  aa  would 
be  expected  with  high  flo*fs.  At  the  very  beginning  and  end  of  the  cheao- 
graph,  deposition  appears. 

The  coaparlsons  between  the  aeasured  and  calculated  chloride  and  conductivity 
are  only  fair.  There  are  several  reasons  for  the  deviation  between  the 
predicted  and  aeaaured  concentrations;  two  are  aajor  reaaona.  First  the 
actual  data  prior  to  the  start  of  the  atora  was  not  used,  the  values  at  the 
atart  of  the  store  were  pre-extended  bock  eo  water  traveling  froa  the 
upetreaa  atatlon  would  reach  the  downstreaa  station  by  the  start  of  the 
store.  TTile  problea  exists  for  aost  of  the  stores  discussed  herein  and  It 
probably  accounts  for  the  alsaatch  at  the  start  of  the  storaa.  Secondly,  the 
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■  Iraatch  occurs  because  Che  unaeasured  Inflow  concent  rat  tons  were  assuaad 
constant  when  In  actuality  they  are  variable.  Other  alnor  contributions  to 
Che  error  could  be  variability  In  the  area-dlscharRe  relationship  or  the 
neplertton  of  dlaperston  of  Che  flowlnjt  water. 

Sandusky  Klver  Wear  Mealco 

Figure  l)a  indicates  a  large  Influence  of  unaeasured  Inflow,  figures  13b~13f 
can  be  exaalned  In  the  aasM  asnner  as  Che  Upper  Sandusky  station.  Figure  13b 
and  Figure  1 3d  Indicate  rasuapenslon,  especially  the  large  resuspenalon  for 
suspended  sedlaent.  Various  concentrations  aust  be  run  to  Indlcste 
unaeasured  Inflow  concentration  Influences.  If  this  Influence  Is  negligible, 
then  these  graphs  realistically  reflect  Che  transport  process  in  Che  systea. 
Generally,  the  predicted  chloride  and  conductivity  (Figures  13e  4  13f)  agree 
with  Che  aeasured  values.  A  variable  unaeasured  lateral  Inflow  concentration 
nay  Inprove  the  conparlaon. 

Sandusky  Klver  Near  Freaont 

Figure  lAa  again  shows  a  significant  contribution  froa  the  unaeasured  basin 
area.  Figures  lAb-lAf  show  a  fairly  good  aatch  for  these  paraaeters  using  a 
single  nonpolnc  chenistry  value.  The  low  values  for  calculated  chloride  and 
conductivity  for  the  22nd  and  23rd  are  due  to  Che  approxlaatlon  of  the  data 
before  the  22nd  using  the  22nd  data.  If  actual  data  ware  used  before  this 
tine,  the  calculated  values  would  be  close  to  the  aeasured  values.  The  total 
phosphorus  Figure  14b  and  s'.ispended  sedlaent  14d  prlaary  Indicate  resuspen- 
slon  with  soae  depoaltlon  especially  at  Che  end  of  Che  atom. 

Llnearltstlon  of  Area-Discharge  Heist lonahlps 

Figure  IS  shows  the  area~dlscharge  relationships  for  the  Bucyrus  (B)  and 
Upper  Sandusky  (US)  stations.  Also  shown  are  two  linearised  versions  of  the 
area-discharge  curve  passing  through  the  origin,  one  labeled  Lj ,  and  the 
other  Lj. 

Figure  16  shows  the  results  obtained  by  routing  Che  hydrograph  at  Bucyrus  to 
Upper  Sandusky  for  a  snail  atom  using  Che  linear  area-discharge  relationship 
labeled  L|.  Since  In  this  case  the  hydrograph  wave  speed  is  the  saae  as  the 
water  speed  and  la  less  chan  It  otherwise  would  be  when  using  the  actual 
area-dlacharge  curve.  It  would  be  expected  that  there  would  be  e  ehlfc  to  the 
right  of  the  hydrograph  properties.  For  this  case  Che  actual  araa-dlscharga 
curves  were  used  to  calculate  N  used  for  Che  above  routing.  Figure  16  can  ba 
coapared  Co  Figure  7  which  uses  the  actual  araa-dlscharga  relatlonahlpa. 

Note  the  aajor  shifts  occur  at  high  flow  where  deviations  of  L|  froa  US  aad  B 
(Figure  13)  are  significant. 

Figures  17a  and  1 7b  show  the  results  for  flow  and  conductivity  for  Upper 
Sandusky  using  linear  area-dlacharge  approx last Ion.  These  results  can  ba 
directly  coapared  to  Figures  12a  and  12f,  respectively.  As  can  ba  aaan  In 
this  case,  there  Isn't  auch  difference  In  the  results  bacauaa  this  llaaarlsa- 
tlon  closely  approxlaatad  the  araa-dlscharga  curves  In  the  flow  region  used. 
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Flgur«»s  I8<i-18c  show  the  results  for  one  day,  February  26,  for  discharge. 
Figure  I8a  Is  the  flow  pattern  for  the  actual  area-discharge  curves,  B  and 
US,  while  Figure  18b  is  for  Lj ,  and  Figure  18c  for  L2.  can  be  seen,  there 
are  slight  differences  between  Figures  18a  and  18b,  while  for  Figure  18c 
there  Is  a  shift  of  alis<>st  half  a  day  In  the  curve  properties  to  tl>e  left. 
There  Is  also  a  change  In  the  'S.»gnltude  of  the  results.  This  shift  and 
changes  In  <«agnltude  can  be  expected  because  the  water  Is  now  forced  to 
travel  faster  than  It  traveled  before  by  using  Lt  Instead  of  8  and  US  In 
Figure  IS;  therefore,  the  water  would  arrive  sooner  downstreaa  and  the  tribu¬ 
tary  and  nonpoint  flows  would  Influence  different  parts  of  the  water  at  Upper 
Sandusky  as  It  traveled  downstream  than  would  have  been  Influenced  before. 

Figures  I9a-I9c  show  the  same  type  of  results  for  conductivity  for  the  same 
day  using  a  nonpoint  concentration  of  *>00  micro  mhos  at  2S°C.  As  can  be 
seen,  th.*  changing  of  the  area-discharge  curves  causes  a  shift  In  the  pre¬ 
dicted  profiles  primarily  because  the  properties  of  the  water  from  different 
sources  Is  shifted  In  time. 


Time  Varying  Inputs 

Figure  20  shows  a  time  varying  min|v>lnt  Inflow  concentration  for  conductivity 
while  Figure  21  shows  the  results.  This  figure  Is  compared  with  Figure  Ilf. 
As  can  he  seen  If  sufficient  Information  Is  known  to  define  a  variable 
lateral  Inflow  concentration,  the  chemistry  may  be  approximated  more  accu¬ 
rately.  The  assumed  profile  In  Figure  20  Is  probably  realistic  since  It 
corresp»>nds  to  the  general  shape  of  the  conductivity  concent  rat  Ion  In  the 
river,  l.e.  Figure  21. 

Figures  22  and  21  show  similar  results  for  total  phosphorus.  Note  the  rela¬ 
tively  Ins  Ignl  f  leant  changes  when  compared  to  Figure  11b.  F.ven  by  Increasing 
the  Inflow  concent  rat  1  on  to  high  values,  there  Is  n<i  significant  changes 
during  the  stt)rm  even  thoiigh  the  nonpoint  Inflow  Is  relatively  large.  Large 
changes  are  soon  only  during  low  flow.  This  Illustrates  that  Inflow  con¬ 
cent  rat  Itnis  cannot  account  for  the  variations  and  therefore,  resuspenslon 
occurs  throughout  this  storm.  At  the  peak  period  of  flow,  the  flow  from 
Upper  Sand\isky  predomi ivttes  and,  because  of  this,  the  cht»mlstry  at  Upper 
Sandusky  predominates  as  set>n  In  the  calculated  peak  during  the  2Sth  as  seen 
In  Figure  lib. 


Instream  Processing  of  Phosphorus 

Figure  24  Illustrates  how  the  calculat lonal  procedure  can  he  used  to 
understand  the  Instream  processing  of  point  source  phosphorus.  This  example 
Is  for  the  Upper  Sandusky  station  which  is  downstream  of  a  aewage  treatment 
outfall,  fbirlng  low  flow,  both  total  and  ortho  phoaphorua  are  loat  from  the 
water  column.  Thirlng  high  flow,  total  phoaphorua  Is  resuspended  while  no 
ortho  phosphorus  Is  removed  from  the  sediments.  This  figure  shows  how  imme¬ 
diately  available  ortho  phosphorus  la  lost  from  the  water  coltssn  and  ts  aub- 
aequently  resuspended  and  delivered  to  the  lake  during  a  atnrm  aa  potentially 
available  pbosphor<is. 
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UPPER  SANDUSKY  TO  MEXICO 
77/02/22 

UNMEASURED  INFLOW  CONCENTRATION 


Ocposttloo  and  Rcauspcnalon  of  («)  Total  Phoaphorua  and  (b)  Orthophoaphata  in  the 
Sanduaky  River  near  Upper  Senduaky  ~  Stora  beginning  7  July  I97g 


A  study  of  different  stonaa  on  the  Sandusky  River  basin  has  Indicated  net 
deposition  to  occur  durlnR  low  steady  floM  and  small  storms  for  aome  reaches 
while  net  resuspenston  occurs  for  large  storms.  The  reach  from  Rucyrua  to 
Upper  Sandusky  generally  shows  net  resuspension  for  all  storms.  Presumedly, 
the  material  that  la  resuspended  Is  deposited  during  low  flows  between 
storms  with  much  material  coming  from  the  sewage  treatment  plant  below 
Bucyrus.  For  the  reach  from  Upper  Sandusky  to  Mexico  net  deposition  occurs 
during  small  storms.  Presumedly,  a  small  dan  above  Mexico  aide  In  the  depo¬ 
sition.  The  reach  from  Mexico  to  Fremont  shows  net  resuspenslon  for  nearly 
all  storms.  Low  steady  flow  causes  net  deposition  In  all  reaches.  An  Impor¬ 
tant  fact  resulting  from  this  study  Is  that  about  half  of  the  total 
phosphorus  entering  Lake  Erie  resulted  from  resuspenslon  of  material  In  the 
main  stream  of  the  river.  Most  of  this  ouiterlal  ««as  Initially  derived  from 
the  land  surface. 


Total  Water  Balance  Check  on  Accuracy 

Table  2  shows  a  comparison  of  the  total  voUae  of  water  for  this  storm  calcu¬ 
lated  by  different  methods  which  Is  a  check  of  the  accuracy  of  the  calcula- 
t tonal  procedure  for  a  particular  storm.  Two  sources  of  numbers  are  depicted 
as  major  colian  headings.  These  are  further  divided  Into  results  depicting 
data  which  exclude  and  Include  the  listed  tributaries  from  the  unmeasured 
lateral  Inflow,  respectively.  The  total  line  Indicates  the  total  of  the 
sources  for  the  particular  reach.  Following  this  Is  the  measured  amount  at 
the  station.  These  two  lines  should  be  In  fairly  close  agreement  for  good 
accuracy.  The  next  line  Indicates  the  stsrtlng  source  for  the  next  reach. 
There  are  two  sources  of  error  attributed  to  these  numbers.  The  flrat,  «d)lch 
Is  small.  Is  due  to  niaierlcal  Integration  of  the  areas  under  the  curves  and 
rounding  to  significant  figures.  The  second  source,  which  could  be  major.  Is 
due  to  downstream  hydrograph  mismatch  as  an  example  shotm  In  Figure  11.  This 
occurs  when  the  tributaries  of  Interest  are  Included  separately  in  the  calcu¬ 
lations  and  their  water  volumes  are  subtracted  from  the  unmeasured  lateral 
Inflow  for  the  reach  with  the  result  being  adjusted  as  not  to  go  negative. 

The  mismatch  volume  Is  used  In  the  calculations!  procedure  with  the  resulting 
percentages  being  applied  to  the  actual  water  volume  to  get  the  numbers  In 
the  colimin  titled  "Calculatlonal  Procedure. **  This  generally  causes  a  minor 
discontinuity  between  stations  (the  total  for  the  reach  and  measured  do  not 
agree).  The  amount  of  error  Is  checked  by  computing  fractional  volumes  id>are 
the  tributaries  are  included  under  the  unmeasured  lateral  Inflow  for  the 
reach.  For  example,  for  the  reach  from  Bucyrus  to  Upper  Sandusky  the  results 
where  the  tributaries  are  not  broken  out  agree  closely.  When  Broken  Sword 
Creek  Is  separated  from  the  unmeasured  inflow  the  total  for  the  reach  Is  760 
as  compared  to  the  measured  Mount  of  733  as  shotm  In  column  1.  Column  3 
shows  the  calculat tonal  procedure  voltsaes  which  are  based  on  the  measured 
flow  of  733.  Because  of  this  there  is  some  error  Introduced  in  the  volumes 
as  shown  when  these  are  compared  to  the  numbers  In  column  A. 

Total  Phosphorus  Balance 

Table  3  shows  a  comprehensible  total  phosphorus  masa  balance  for  the  Sanduaky 
River  for  the  storms  of  77/02/22  which  lasted  for  10  days  as  derived  by 
finding  the  area  under  the  total  phosphorus  fraction  graphs  such  aa  the  one 
In  Figure  12b  derived  by  fraction  calculation.  The  budget  must  be  viewed  In 
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Table  2  -  Coaparlaon  of  the  Total  Voluae  of  Water 
Stora  of  77/02/22 

(All  voliaaea  are  In  ■llllona  of  cubic  feet) 


Sourcea  -  Froa 

Rucyrua  to  Freaont 

Froa  Input  Data 

Calculatlonal 

Procedure 

Bucyriia 

350* 

350* 

321 

352 

Broken  Sword  Creek 

267 

- 

268 

- 

Unaeaaured 

143 

383 

144 

382 

Total  (Bucyrua  to  Upper  Sandusky) 

760 

733 

733 

733 

Upper  Sandusky* 

733 

733 

733 

733 

Upper  Sandusky 

730* 

730* 

740 

740 

Tyaorhtee  Creek 

350 

- 

340 

- 

Unaeaaured 

1,200 

1,550 

1,190 

1,540 

Total  (Upper  Sandusky  to  Mexico) 

2,280 

2,280 

2,270 

t 

2,280 

Mexico* 

2,270 

2,270 

2,270 

2,270 

Mexico 

2,270* 

2,270* 

2,110 

2,230 

Honey  Creek 

450 

- 

450 

- 

Wolf  Crek 

1,060 

- 

1,040 

Unaeaaured 

1,650 

3,030 

1,690 

3,060 

Total  (Mexico  to  Freaont) 

5,430 

5,300 

5,290 

;  5,290 
: 

:  5,290 

Freaont* 

5,290 

5,290 

t 

5,290 

*H^MMured  field  data  and  uaed  aa  Input  to  the  calculatlonal  procedure 
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Table  3  -  Kxaaple  of  a  Detailed  Total  Phosphorus  Hass  Balance 

Stora  of  77/02/22 


Sources  -  Fro«  Bucyrus 

Flow 

(allllons) 

Total  Phosphorus 
(aetrlc  tons) 
(unmeasured  Inflow 
concentration  0.2  bk/I) 

Bucy rus 

321 

6.22 

Broken  St#ord  Creek 

268 

2.21 

Unaeasured 

144 

0.81 

Net  Resuspenslon/Deposlt  Ion 
(Bucyrus  to  Upper  Sandusky) 

- 

12.63 

Total  (Bucyrus  to  Upper  5>andusky) 

733 

21.87 

Upper  Sandusky 

740 

22.09 

Tyaochtee  Creek 

340 

3.03 

Unaeasiired 

1,190 

6.75 

Net  Resuspenslon/Deposlt Ion 

. 

(Upper  5tanduaky  to  Mexico) 

- 

24.01 

Total  (Upper  Sandusky  to  Mexico) 

2,270 

55.88 

Mexico 

2,110 

54.  78 

Honey  Creek 

450 

14.14 

Wolf  Creek 

1,040 

4.25 

Unaeasured 

1,690 

9.56 

Net  Resuspenslon/Deposlt  ion 
(Mexico  to  Freaont) 

- 

17.24 

Total  (Mexico  to  freaont) 

(for  Freaont) 

5,290 

99.97 

\ 


Table  4  -  Exaaplc  of  a  General  Total  Phoaphorua  Maaa  Balance 
Stora  of  22  February  1977 


th»  conCext  of  the  above  dlscusalon  of  the  discontinuities.  The  sources 
traced  In  the  Model  can  be  few  or  many  depending  upon  the  details  desired  and 
the  availability  of  the  data.  The  corresponding  budgets  can  be  detailed  or 
general  as  desired.  Table  4  shows  a  general  total  phosphorus  budget  where 
the  calculation  discontinuities  between  reaches  In  Table  3  have  been  smoothed 
out.  Similar  calculations  could  be  performed  for  the  parameters,  e.g. 
suspended  sediment  or  chlorides,  to  get  a  suspended  sediment  or  chloride 
stream  budget. 

Management  Implications 

By  applying  the  ca Iciilat Iona  1  procedure  to  various  storms  an  engineer  or 
planner  could  develop  land  or  stream  management  practices.  The  calculations 
may  be  used  to  Indicate  a  problem  subbasln  area  which  appears  for  all  storms. 
In  the  cane  of  the  5>anduaky  River  as  stated  before,  about  half  the  swasured 
total  phosphorus  results  from  resunpension  of  bottom  materials  which  orlgl~ 
nated  from  land  runoff.  Presumedly  such  material  Is  carried  %rlth  the  sedi¬ 
ment  eroded  from  the  land.  Improvements  In  erosion  control  and  changes  In 
land  use  will  change  the  amount  of  material  that  shows  up  In  the  stream. 
Immediate  changes  In  water  quality  may  be  masked  by  the  randomness  of  the 
phenomenon  that  we  are  dealing  with.  True  changes  may  not  be  evident  for  a 
couple  of  years.  Point  source  Inputs  tend  to  deposit  during  low  flow  pro¬ 
viding  a  source  of  material  for  resuspension  during  high  flows.  In  such 
cases  available  phosphorus  Is  "hanked,"  transported  by  resuspension,  and 
delivered  as  potentially  available  total  phosphorus.  Besides  controlling 
sediment  and  phosphorus  at  its  source  on  the  land  It  may  be  possible  to 
control  them  in  the  stream  by  structural  means.  By  using  retention  struc¬ 
tures  material  could  be  forced  to  settle  during  low  flows  and  retained  there 
during  high  flows,  thus  limiting  the  amount  of  Input  for  resuspension. 

CONCLUSION 

The  calculation  of  conservative  substances  (chloride  and  conductivity) 
generally  Indicate  the  calculations!  procedure  developed  herein  works  %(ell  In 
spite  of  the  various  approximations  used.  Using  variable  unmeasured  Inflow 
concentrations  during  storms  for  conservative  substances  give  good  matches 
for  the  measured  snd  predicted  values.  Suspended  sediment  and  total 
phosphorus  usually  show  resuspension  for  rising  flows  and  deposition  for 
declining  flows.  This  results  in  a  net  overall  resuspension  for  the  storm. 
This  withdrawal  Is  balanced  by  a  net  deposition  during  low  flow  periods. 
Resuspension  of  total  phosphorus  during  storm  events  la  shorn  downstream  from 
a  phosphorus  point  source.  Total  phosphorus  and  other  phosphorus  are  depo¬ 
sited  from  the  stream  during  low  flow  conditions  and  only  total  phosphorus  Is 
removed  from  the  sediments. 

SUMMARY 


A  calculations!  technique  has  been  presented  for  tracing  the  transport  of 
nutrients,  suspended  solids,  and  other  chemical  substances  in  a  stream.  At  a 
given  point  in  a  stream,  these  auiterlsls  are  derived  from  various  upstream 
sources.  This  procedure  permits  the  calculation  of  the  fractions  of  material 


derived  froa  these  sources.  Also  Che  calculation  of  the  aaounC  of  resuspen- 
alon  and  deposition  of  total  phosphorus  and  suspended  solids  In  various 
streaa  reaches  Is  possible. 

Althoiifih  the  calculation  of  conservative  substances  Indicates  Che  aethod 
works  fairly  well  In  spite  of  Che  various  approxlaat Ions  used,  It  aust  be 
reaeabered  Chat  this  procedure  was  derived  for  the  prlaary  purpose  of  Indi¬ 
cating  what  Is  happening  at  the  downstreaa  station  of  a  reach  during  a  stora 
event  and  not  for  Che  prediction  of  Che  chealstry  along  the  reach. 

It  Is  also  possible  to  develop  total  nutrient  auss  balances  Indicating 
contributions  froa  various  areas  of  the  watershed. 

This  procedure  was  developed  with  Che  objective  of  providing  Insight  as  to 
what  Is  happening  to  various  sources  of  Input  aaterlals  with  Che  least  aaount 
of  Input  data  required.  It  la  hoped  this  technique  will  provide  researchers 
with  a  aethod  of  evaluating  InsCreaa  processing  of  various  aaterlals  froa 
various  watershed  areas. 

Hore  details  can  be  obtained  by  referencing  references  5-12  and  15. 
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